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Abstract 

We study self-contact phenomena in elastic rods that are constrained 
to lie on a cylinder. By choosing a particular set of variables to describe 
the rod centerline the variational setting is made particularly simple: the 
strain energy is a second-order functional of a single scalar variable, and the 
self-contact constraint is written as an integral inequality. 

Using techniques from ode theory (comparison principles) and varia- 
tional calculus (cut-and-paste arguments) we fully characterize the struc- 
ture of constrained minimizers. An important auxiliary result states that the 
set of self-contact points is continuous, a result that contrasts with known 
examples from contact problems in free rods. 
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1. Introduction 



The study of self-contact in elastic rods has seen some remarkable progress 
over the last ten years, with highlights such as the numerical work of To- 
bias, Coleman, and Swigon [22,6,5], the introduction of global curvature 
by Gonzalez and co-workers [10], and the derivation of the Euler-Lagrange 
equations for energy minimization by Schuricht and Von der Mosel [19]. 
Parallel advances have been made on the highly related ideal knots and 
Gchring links, where ropelength is minimized instead of elastic energy [4, 
18,3]. 
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Despite this progress important questions remain open. We are still far 
from understanding analytically the solutions of the Euler-Lagrange equa- 
tions for general contact situations. Even if we limit ourselves to global 
minimizcrs of an appropriate energy functional, we can prove little about 
the form of solutions as soon as contact is taken into account. 

For instance, a long-standing conjecture for closed elastic rods is that 
in the limit of long rods under constant twist the global energy minimizer 
should be a ply (double helix) with a loop on each end. If a structure of this 
type is assumed, then the limiting pitch angle can be determined [21]; but 
the difficult part actually consists in showing that global minimizers have 
this structure. Incidentally, since local minimizers of different type have 
been found numerically [6,5], the restriction to global minimizers appears 
to be essential. 

This example is typical for the current state of understanding: if assump- 
tions are made on the set of contacts, then characterizations are possible [8, 
20, 21, 12], but for unrestrained geometry little is known rigorously. It shows 
how our lack of understanding of energy minimizers is intimately linked to 
the lack of knowledge about structure of the contact set. Examples show 
that this structure can be non-trivial: for instance, non-contiguous contact 
appears at the end of a ply in an elastic rod [5] . 

In this paper we study a problem of self-contact of an elastic rod in 
which the rod has reduced freedom of movement: the centerline of the rod 
is constrained to lie on the surface of a cylinder (Figure 1). In contrast to the 
full three-dimensional problem referred to above, the reduced dimensionality 
of this problem enables us to give a near-complete characterization of global 
minimizers, without making any a priori assumptions on the structure of the 
contact set. Notwithstanding this, determining the structure of the contact 
set is a central element of this paper. 

contact point 



x 




Fig. 1. The centerline of a rod on a cylinder is described using cylindrical coordi- 
nates: the independent variable x is the tangential coordinate, and the position of 
the centerline is given by the function ((x) measuring distance along the cylinder 
axis. 

We transform the classical Cosserat model of an elastic, unshearable rod 
of circular cross-section into a more convenient form. The functional that 
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is to be minimized (representing stored energy and work done by the end 
moment) is 



F(u)= f [o(u)u /2 + &(«)]. 
Jo 



where 



°(") = 7Z9 ^ , ..^ R /9 and 



4tt 2 (1 + u 2 )5/ 2 w r 2 (l + w 2 ) 3 / 2 Br VTT 
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(1) 

Here r is the radius of the cylinder, M the moment applied to the end of 
the rod, and B is the bending coefficient of the rod. The centerline of the 
rod is characterized by £(x), which measures distance along the cylinder 
axis as a function of a tangential independent variable x. The unknown in 
this minimization problem is the derivative u(x) — (,'{x), which may by 
thought of as the cotangent of the angle between the centerline tangent 
and the cylinder axis; u is zero when the rod curls around the cylinder 
orthogonal to the axis, and u = ±oo when the rod is parallel to the axis. 
This transformation is detailed in Section 3. 

The most interesting part of the variational problem is the transformed 
contact condition (condition of non-self-penetration) . In this paper we take 
the thickness of the rod to be zero; then the non-self-penetration condition 
is 

fX+l 

u>0 for all 0<a:<T-l, (2) 



/ 

J X 



where the interval [x, x + 1] corresponds to one full turn around the cylin- 
der; this condition formalizes the intuitive idea that non-self-penetration is 
equivalent to the condition 'that the rod remain on the same side of itself. 
This condition on u makes the variational problem a non-local obstacle prob- 
lem. Non-zero thickness requires a contact condition that is substantially 
more involved than (2); we comment on this situation in Section 3.5. 

Both the background in rod theory and the independent mathematical 
context of this minimization problem raise questions about the solutions: 

1. Do solutions exist? 

2. What is the minimal, and what is the maximal regularity of minimizers? 

3. When is there contact, i.e., when is the contact set 

w c := jxe [0,T-1] : J u = oJ (3) 

non-empty? 

4. Given that co c ^ 0, what is the structure of w c ? Is the contact simply 
contained in a single interval, or is the structure more intricate, as in 
the examples of contact-skip-contact at the end of a ply [6] and in a 
(ropelength minimizing) clasp [3]? 

5. What form do the contact forces take? 
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6. Does the solution inherit the symmetry of the formulation? This is the 
case for a symmetric rod on a cylinder without contact condition [11], 
but need not be true when taking contact into account. 

In the rest of this paper we address these questions. 



2. Results 

The first main result of this investigation (Theorem 1) shows that the 
contact condition (2) is essential: without this condition the centerline of 
a rod will intersect itself. A little experiment with some string wrapped 
around a pencil will convince the reader that this is the case. We also prove 
the regularity result that a constrained minimizer u is of class W 2 ' 00 , and 
we derive the Euler-Lagrange equation 

N(u){x) :=-2a(u(x))u"{x)-a'{u{x))u' 2 (x)+b'(u(x)) = f /, (4) 

Jx-1 

where the Lagrange multiplier / is a non-negative Radon measure with 
support contained in the contact set uu c (Theorem 2). 

From stationarity alone, which is the basis of Theorem 2, the character- 
ization of / as a positive Radon measure appears to be optimal; no further 
information can be extracted. In Section 7 we use two different additional 
assumptions to further characterize the contact set and subsequently the 
measure /. In both cases we obtain the important result that the contact 
set is a (possibly empty) interval and that the measure / is a sum of Dirac 
delta functions, as represented schematically in Figure 2. The weighting of 
the delta functions is shown in the middle of Figure 2: there is a linear 
decrease or increase in weight from one side of the contact set to the other 
(Theorem 3). Since / may be interpreted as the contact force, we deduce 
that 

— The contact force is concentrated in at most two tangential positions X\ 
and X2-, and in integer translates of £1,2! 

— The magnitude of the contact force is maximal at the contact point 
where the rod lifts off, and decreases linearly with each turn. Figure 3 
graphically illustrates this behaviour. 

The decrease in contact force with each turn can be understood in the 
following way. The difference between the contact forces on either side of 
the rod creates a resulting force exerted on the rod, and the two resultant 
forces that act at x\^ mod 1 point in opposite directions. If we imagine a 
single, closed ring with two forces acting on it in this way, the two forces 
create a moment that will bend the ring. This also happens with the coil of 
the current problem, as is demonstrated by the small but definite oscillations 
in the numerical solutions calculated in Section 9. 
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Fig. 2. The function g(x) = j x _ 1 f is piecewise constant (top); the jumps cor- 
respond to Dirac delta functions in / (middle). Note that the support of g is the 
set lo c + [0, 1] by the definition of g. The solution u corresponding to / and g is 
shown at the bottom. 






Fig. 3. A typical rod configuration (left; front and back views) that minimizes 
energy and satisfies the contact condition. On the right the bars indicate the con- 
tact forces corresponding to the arrows in Figure 2. (The analysis of this paper 
assumes zero rod thickness — in this picture the rod has been fattened for presen- 
tation purposes.) 



As mentioned above, the crucial result that the contact set is connected 
requires additional assumptions. If we step back from this rod-on-cylinder 



6 



Van der Heijden, Peletier, and Planque 



model, and allow a and b to be general given functions, then for a large class 
of such functions the nonlinear operator on the left-hand side of (4) N(u) 
satisfies a version of the comparison principle, 

Nui > Nu 2 =^> u\ > U2, 

(see Definition 2 for the precise statement). For such functions a and b, any 
stationary point has a connected contact set (Theorem 4). The argument 
is based on the observation that non-contact in some interval (a, (3) implies 
that / = on {a, (3) and therefore that the right-hand side of (4), 

9{x) = f f, (5) 

Jx-l 

is non-increasing on (a, (3) and non-decreasing on (a — 1, (3 — 1). 

Importantly, however, the functions a and b given in (1) are such that 
the associated operator mostly fails to satisfy this comparison principle. 
We therefore also take a different approach, in which we obtain the same 
result by only considering global minimizers. Here we use a different argu- 
ment, based on constructing other minimizers by cutting and pasting; the 
combined condition of minimization and non-contact in an interval (a, (3) 
implies the existence of additional regions of non-contact outside of the in- 
terval (a, (3), implying that the right-hand side of (4) is constant on (a, (3). 
From this the result follows (Theorem 5). 

In both cases, the fact that the contact set is an interval implies that the 
boundary of the contact set is 'free' — the measure / is zero on an additional 
interval of length one extending on both ends of u> c . This implies that the 
right-hand side g is increasing and decreasing at the same time — except 
at points that lie at integer distance from the two boundary points. This 
imposes the specific structure on g and / that is shown in Figure 2. 

The issue of symmetry of minimizers is a subtle one, which again de- 
pends on the presence or absence of a comparison principle. The compari- 
son principle simplifies the structure of solutions: all stationary points are 
symmetric (up to an unimportant condition on b). Without a comparison 
principle, and more precisely when minimization of F favours oscillation, 
this is no longer true, and even stationary points that are global minimizers 
may be asymmetric (Section 8). 

Using the characterization of / and g derived earlier we use two nu- 
merical methods to investigate constrained minimizers (Section 9); one is a 
method of direct solution, using a boundary-value solver, and the other a 
continuation method. A typical solution is shown in Figure 3. 

The simple structure of the functional and the contact condition suggest 
that the methods and results of this paper might be applicable to other 
systems than this particular rod-on-cylinder model. We therefore state and 
prove our results for general functions a and b. The main requirements 
are that a and b are smooth and that a is positive; other conditions are 
mentioned in the text below. 
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3. Problem setting: derivation of the rod-on-cylinder model 

3.1. Kinematics 

Consider an clastic rod of circular cross-section that is constrained to lie 
on a cylinder, and which is subject to a force T and a moment M at the ends. 
We assume that at the rod ends, T, and M are maintained parallel to both 
the axis of the cylinder and the axis of the rod, but that the loading device 
leaves the rod ends free to rotate around the circumference of the cylinder; 
the ends of the rod therefore need not be coaxial. The rod is naturally 
straight and inextensible, and material cross-sections are assumed to remain 
orthogonal to the centerline. We will derive a minimization problem for rods 
of length 21 and later take the limit I — > oo. 

In the Cosserat rod theory [1, Ch. VIII] the configuration of this rod is 
characterized by a right-handed orthogonal rod-centered coordinate frame 
of directors, {di,d 2 ,ds}, each a function of the arc-length parameter s. 
The director d 3 is assumed parallel to the centerline tangent, and by the 
assumption of inextensibility the centerline curve r satisfies 

r = d 3 , 

where the dot denotes differentiation with respect to arc length. The strain 
of the rod is characterized by the vector- valued function u given by 

dk = u x dk, k = 1, 2, 3. 

When decomposed as u — K\d\ + n 2 d 2 + rd 3l the components may be 
recognized as the two curvatures n\^ 2 and the twist r. 

We choose a fixed frame of reference {ei, e 2 , e%}, where e$ is parallel to 
the cylinder axis, and we relate the frame {di,d 2 ,d 3 } to this frame by a 
particular choice of Euler angles {9, ip, (j)} [14, 11]. In this parametrization 9 
is the angle between d^ and e% (or between the centerline and the cylinder 
axis), V characterizes the rotation around the cylinder axis, and is a 
partial measure of the rotation between cross-sections. The condition that 
the centerline of the rod lie on the surface of a cylinder of radius r translates 
into the kinematic condition 

tp = — sin 9, (6) 
r 

where the dot denotes differentiation with respect to the arclength coor- 
dinate s. Note that it is natural not to restrict tp to an interval of length 
2tt. In terms of the remaining degrees of freedom {9, (/>} the curvatures and 
twist are given by 

K\ = 9 sin 4> — — sin 2 9 cos <f), 
r 

1 „ 

k 2 = 9 cos (f> H — sin 9 sin 0, 
r 

t = <f> + — sin 9 cos 9. 

r 
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3.2. Energy, work, and a variational problem 



For a given rod the strain energy is given by [14], 




4 * + ?// 



Here B and C are the bending and torsional stiffnesses respectively. To 
determine the work done by the tension and moment at the ends of the 
rod we need to characterize the generalized displacements associated with 
these generalized forces. For the tension T the associated displacement is 
the shortening S, 



The generalized displacement associated with the moment M is the end 
rotation, which is well-defined by the assumption of constant end tangents. 
It is common to identify the end rotation with a link-like functional 



As demonstrated in [13], this identification is correct in an open set around 
the undeformed configuration 9 = 0, but loses validity when |0| takes values 
larger than it. Although nothing we have seen suggests that in an energy- 
minimizing situation 9 would take values outside of the admissible interval 
(— 7r, 7r), we have no rigorous argument to guarantee that 9 remains inside 
that interval, and therefore we are forced to assume this. In terms of the 
variables 9 and t this functional then takes the form 



Here we assume rigid loading in shortening and dead loading in twist, 
i.e. we prescribe the shortening S and the moment M , which implies that 
the tension T and the end rotation L arc unknown and to be determined 
as part of the solution. This loading condition leads to the minimization 
problem 



for given a > 0. The tension T has a natural interpretation as a Lagrange 
multiplier associated with the constraint of S. 

We can simplify this minimization problem by first minimizing with 
respect to r for fixed 9, from which we find r = M/C; re-insertion yields 
the final minimization problem 




L = 





min {E{9, r) - ML(9, r) : S(9) = a} 



min {F(9) : S(9) = a} 



(7) 
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with 

F(0) = | Jo 2 + 7^ J^9-^- /sin0(l-cos0). (8) 

We are interested in localized forms of deformation, in which the defor- 
mation is concentrated on a small part of the rod and in which boundary 
effects are to be avoided, and therefore we take an infinitely long rod and 
consider 9, F, and S to be defined on the whole real line and assume 9 — > 
as \s\ — > oo. 



3. 3. Behaviour of minimizers 

The Euler-Lagrange equations associated with the minimization prob- 
lem (7) can be written as a Hamiltonian system with one degree of freedom, 

l -9 2 + V{6) = H, (9) 

for a particular V. In this system two independent parameters remain, which 
may be interpreted as a scaled cylinder radius f = rM/B and a combined 
loading parameter m = Mj \J BT. 

Solutions of the original minimization problem are orbits of this Hamil- 
tonian system that are homoclinic to zero, and such orbits have been studied 
in detail in [11]. Among the findings are 

1. For all values of f ranges of m exist with orbits that are homoclinic to 
the origin; 

2. At some parameter points these homoclinic orbits 'collide' with saddle 
points. The saddle points correspond to helical solutions (constant angle 
9) and close to these collisions the homoclinic orbit has a large region of 
near-constant angle 9. 

In Figure 4 a bifurcation diagram is shown with two such collisions, one at 
a forward helix (0 < 9 < tt/2, at m — m C2 ) and one at a backward helix 
(ir/2 < 9 < it, at m — m Cl ). 

In [11] the question of stability of these solutions, both local and global, 
was left untouched. If we interpret the combined load parameter m as 
a (reciprocal) tension T (with the moment M fixed) then the nature of 
the bifurcation diagram in Figure 4, involving as it docs the mechanically 
conjugate variables S and T, suggests that in each peak the right curve is 
locally stable [15]. With two peaks occurring however, this does not allow 
us to predict where the globally stable solution is located. 

In this paper we focus on global energy minimization. Corollary 1 be- 
low states that for sufficiently large shortening, and when contact effects 
are neglected, global energy minimizers always intersect themselves. It is 
this result that forms the main motivation of the analysis of this paper: 
since energy minimization without appropriate penalization leads to self- 
intersection, the non-self-intcrsection condition is necessary for physically 
acceptable solutions. 



10 



Van der Heijden, Peletier, and Planque 




(d) (e) (f) (g) 

Fig. 4. A load-displacement diagram showing shortening d — SM/B of station- 
ary points as a function of the (combined) load m (from [11]). Contact effects are 
not taken into account. The peaks divide this diagram into three sections. The 
solutions in the middle section intersect themselves, whilst the solutions on the 
right do not. The section on the left consists of heteroclinic connections between 
9 = and 6 = 2-k which are not considered here. For sufficiently large shortening, 
the rod configuration that has lowest energy is on the self-penetrating branch, as 
shown by Corollary 1. 

3.4- Translation to (u,ip)- coordinates 

To study the case in which self-contact is taken into account, it is neces- 
sary to properly restrict the class of admissible functions in the minimisation 
problem (7) . In three dimensions a variety of different descriptions of self- 
contact exists for rods of finite thickness, each with subtle advantages and 
disadvantages (see e.g. the introduction of [10]). For a rod on a cylinder 
the situation is simpler, since the freedom of movement is essentially two- 
dimensional — similar to that of a curve in a plane. We focus on rods of 
zero thickness, and implement non-self-penetration as non-self-intersection 
of the centerline. In terms of the unknown #(•) as introduced above, this 
condition can be written as 

z(si) — z{s 2 ) 7^ for all S\ ^ s 2 with tp(si) — -0(s 2 ) = mod 2ir, 

(10) 

where we have used the previous equation (6) for ip and the axial coordinate 

z: 

if) = - sin0, z = cos 9. 
r 

We now make the assumption that z can be written as a function of ip, or, 
equivalently, that ip is monotonic along the rod. This assumption is satisfied 
for solutions of the problem without contact having 9 < w, as given by 
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equation (9). If we include a contact condition of the form (10), then we are 
unable to prove that ip is monotonic, and in fact it is conceivable that this 
monotonicity is only valid for global energy minimizers. 



Under the assumption that z can be written as a function of ip, we intro- 
duce a dimensionless axial coordinate £ = z/r, and write ' for differentiation 
with respect to tp. The functional F in (8) then transforms to 



no 



Br 



C" 2 



(i + c¥ 



B 
Yr 



(i + c¥ 



M 



yj + C 2 - C 



with shortening 



s(0 



i + C' 2 - c 



Here [0,T], the domain of definition of ip, is a priori unknown, since the 
ends of the rod are free to move around the cylinder. 

In these variables non-self-intersection is easily characterized. Since ip is 
monotonic, let us assume it to be increasing (this amounts to an assumption 
on the sign of the applied moment M). Admissible functions arc defined by 
the following condition: 



V?/> e [0, T - 2ir] : Q(ip + 2n) - ((ip) > 0. (11) 

Note that it is only necessary to rule out self-intersection after a single 
turn; if contact exists after multiple turns, contact also exists (potentially 
elsewhere) after a single turn. 

The contact condition (11) is the novel part in this variational prob- 
lem. In this paper we focus on the effect that this condition has on the 
minimization problem, and therefore simplify by 

— fixing the domain size T, and accordingly removing the shortening con- 
straint; 

— replacing the mechanically correct boundary conditons (' = oo by a 
more convenient condition £' = 1. 

In terms of the new variables x = ip/2n and u(x) — C'(V') we recover the 
problem of the introduction. 

These boundary conditions can be described as follows. By prescribing 
£' = u = 1 at the ends of the rod we fix the angle between the rod and 
the centerline to 7r/4. By removing the shortening constraint we allow the 
ends of the rod to move freely in the axial direction; in contrast, the fixing 
of T prevents the rod ends from moving tangcntially. We believe that these 
changes have little effect on that part of the rod that is implicated in the 
contact problem; but this is a topic of current research. 
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3.5. Zero thickness 

The assumption of zero rod thickness can not be relaxed without intro- 
ducing important changes in the formulation (see Figure 5). At thickness e, 
the distance in the ^-direction between two parallel consecutive centerlincs 
in contact is e/sinfl, where 9 is the angle between the centerlines and the 
cylinder axis. Therefore non-zero thickness can not be introduced by simply 
replacing the right-hand side in (11) by e; the angle of the centerlines is to 
be taken into account, implying that the right-hand side of (11) will depend 
on C- 

To make matters worse, when the centerlines are not parallel, i.e. when 
u = (' is not constant, the minimal-distance connection between two consec- 
utive turns depends on values of £' nearby (see [16] for a thorough treatment 
of the geometry of this issue) ; it is not clear whether for the present case 
of a rod on a cylinder any simpler impenetrability condition can be found 
than the well-known global curvature condition [10]. 




Fig. 5. Two configurations of a rod of thickeness e. This illustrates that for rods 
with positive thickness one cannot simply replace the contact condition Bu > 
by Bu > e; a more involved condition is necessary. 



4. Existence and the contact condition 

In this section we state precisely the problem under discussion and show 
that minimizers exist. We also study the minimization problem without the 
contact constraint, and show that minimizers will intersect themselves. 

Let U = 1 + X, where X = H£(0,T), and Y = C([0,T - 1]). Let the 
functional F : U — > R be defined as in the Introduction, 

f{u)= / [o(t*y 2 + &(«)]. 

Jo 

and introduce the constraint operator B : U — > Y given by 

rx+l 

B(u)(x) = u. 

J x 
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With the set of admissible functions given as 

K :={ueU : B{u){x) > Vx e [0, T - 1]} 

the central problem is 

Problem (A): Find a function u* <G U such that 

F{u) = min{F(u) : u e K}. 

We first prove existence of minimizcrs for Problem (A). 

Lemma 1. Let T > 0. Assume that a(u) > ao > 0, and that b(u) is Lips- 
chitz continuous. Then there exists u* e K such that 

F(u*) = min{F(u) : ueK}. 

Proof. Let {u n } C K C U be a minimizing sequence. We first prove that 
/ b(u n ) is bounded from below. 

Since minimization of F is equivalent to minimization of F — T6(0), we 
can assume without loss of generality that 6(0) = 0. Using the Lipschitz 
continuity of b and the Poincare inequality we have 

\\b(u n )\\ L1 <c\\u n \\ L1 < Cl (T+K-l|| L1 ) < C (1 + ||<|| L2 ). 

Here and below c is a possibly changing constant that does not depend on 
n. Then 

J b{u n ) > 

> 
> 
> 

Hence 



for a suitable constant D 

Together with the boundedness of F(u n ), (12) implies that u n is bounded 
in X. Hence {u n } contains a subsequence {u„ m } that converges weakly in 
X to a limit u* . Since F is lower semicontinuous with respect to weak 
convergence, 

F(u*)< liminf F(u nm ), 

m—>OQ 

implying that W* is a minimizer. 



-c(i + KM 
_c ( 1 + i(/ aKR2 ) 2 ) 

-c(l + ^(F{u n ) - J b{u n )) h ^j 
- c (l + l( c -/ &K ))'). 

Jb(u n )>-D (12) 
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As we mentioned in Section 3, if contact is not taken into account — if 
F is minimized in U rather than in the smaller set K — condition. In the 
theorem below we actually prove a stronger statement. We write Ft instead 
of F to indicate explicitly the dependence on the interval [0, T]. 

Theorem 1 (Minimization without contact condition). Assume that 
a and b are of class C 1 , and that a is strictly positive. Assume that some 
u < 1 exists such that 

-oo < inf 6 < inf b. (13) 

R u>u 

For each T > 0, let ut be a minimizer corresponding to the minimization 
problem on domain [0, T], 

min{F T (M) : u e U}. (14) 

Then there exists c > independent of T such that 

\{xe [0,T] :u T {x) >u}\ <c(l + Vf). 

The function b given in (1) achieves its minimum at u = — oo, regardless 
of the value of Mr/B; therefore it satisfies the condition (13) for every 
u < 1. 

Corollary 1 (Minimizers violate the contact condition). In addition 
to the conditions of Theorem 1, assume that u < 0. IfTis sufficiently large, 
then B(u T ){x) < for some x e [0,T - 1]. 

Proof (Proof of Theorem 1). We first use a standard argument to give 
an upper bound on the energy Ft(ut)- Choose a T-dependent constant 
u T < u such that 

b(u T ) < inf b(u) and < b(u T ) - inf b < T" 1/2 . 

u>u R 

For any T G M + we construct a new continuous symmetric function ut G U 
such that ut = ILt on — 1]> an d Ft(ut) < C + Tb(u), where C does 
not depend on T. Since ut minimizes Ft, it also follows that 

F T {u T ) < F T {u T ) < C + Tb{u T ). (15) 

Among other things this inequality implies that for large T a minimizer ut 
can not be the constant function 1. 

The Eulcr-Lagrangc equation associated with this minimization problem 

is 

-2a(u)u" - a'(u)u' 2 + b'(u) = 0, (16) 
which can also be written als a one-degree-of-freedom Hamiltonian system 

-a(u)u' 2 + b{u) = H. (17) 

It follows that for any minimizer u, 
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1. b(u(x)) = H at any stationary point x of u; 

2. b(u(x)) > H for all x G [0,T]; 

3. 6(1) > H. 

The third statement follows from noting that if 6(1) = H then u = 1 would 
be the unique solution of (17). 

We now show that any minimizer u is bi-monotonic, i.e. increasing or 
decreasing away from a minimum or maximum. Suppose instead that u 
has two internal stationary points, a minimum at x± and a maximum at X2, 
assume for defmiteness that < x\ < x 2 < T. Note that u{x\) < 1 < u{x 2 ) 1 
since the solution of the Hamiltonian system is a periodic orbit oscillating 
between the values u{x\) and u(x2); the inequality u{x\) < 1 < u{x 2 ) 
follows from the boundary condition. Now pick a point x\2 G {x\,x 2 ) such 
that u{x\2) = 1- 

Construct a new function 

!u(x) < x < xi 

u(xi) xi < x < xi + T — x\2 

u(x-T + xi 2 ) xi+T -xi 2 <x<T 

Then 

F T (u)= { \a{u)u' 2 + b(u)}+ f 1+ " b{u{xi))+ f ^ [a(u)u /2 + b(u)} 

JO J xi J xi 

fXl2 f>Xi+T—Xi2. 

= / [a(u)u' 2 + b(u)]+ / b(u(xi)) 

JO J xi 

= f ^ [a{u)u' 2 + b{u)] + H{T - xi 2 ) 
Jo 

< F T {u). 

Therefore the assumption of two stationary points is contradicted. Note that 
by (17) the solution also is symmetric in [0,T]. 

We now return to the sequence of functions ut- Setting A = {x e [0, T] : 
ut(x) > u} we have 

C + Tb{u T ) > F T {u T ) 

>\A\ inf 6( U ) + (T-L4|)inf& 

u>u R 



U| ( inf b(u)-mib) +Tmib, 

\u>u R / R 



so that 



Ml ( inf b(u) -Mb) <C + T(b(u T ) - inf 6) < c (1 + VT). 

\u>u R / R 

This concludes the proof. 

Remark 1. By a very similar argument one may show the following state- 
ment: if minn 6 is uniquely achieved at some u e R, then 



\\ut - u \\l^(Vt,t-Vt) 



as T — > oo. 
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5. The Euler-Lagrange equation 

We characterize the duality (X, X') by identifying the smooth functions 
on [0, T] with a dense subset of X' via the duality pairing 

x>(Z,x)x = / 
Jo 

Similarly we identify Y' with the space of Radon measures RM([0, T — 1]) 
via the same duality pairing, defined for smooth functions, 

y'(v,v}y = / w- 
Jo 

Where necessary, we extend Radon measures in Y' by zero outside of their 
domain [0,T- 1]. 

Theorem 2. Assume that a and b are globally Lipschitz continuous, and 
that a > a > 0. Let u g U be a solution of Problem (A). Then u g 
TT 2,oo (0,T) and there exists a Radon measure f GY' such that 

-2a(u(x))u" (x) - a'(u(x))u' 2 (x) + b'(u(x)) = f f(s) ds (18) 

Jx-l 

for almost every x g (0,T). Moreover f > and supp / C u> c . 

Definition 1. A function u g [/ is called a stationary point if it there exists 
a Radon measure / G V, with / > and supp / C u c , such that (18) is 
satisfied. 

In the rest of the paper we will often drop the arguments in (18) and write 
-2a(u)«" - a'(u)u' 2 + b'(u) = f f. 

Jx-l 

The proof of Theorem 2 follows along the lines of [2] . We fix the function 
u, with contact set uj c defined in (3), and introduce the cone of admissible 
perturbations V, 

V := {ve X : 3{e n } neN C R + ,e„ -> such that B(u+e n v) > OVn g N}. 

Lemma 2. Let u be a minimizer. Then F'(u) ■ v > for all v g V. 

Proof. For any v £ V, the fact u is a minimizer implies that 

F(w + e n v) - F(u) > for all n g N. 

The conditions on a and 6 imply that F is Frechet diffcrcntiable in it (this 
follows from the conditions on a and inspection of (20) below), so that 

< F{u + e n v) - F(u) = e n F'{u) ■ v + o(e n ||t>||x), 

from which it follows that F'(u) ■ v > 0. Now, given any v g V, take a 
sequence v m C V that converges to w in X. Since F'(u) : X — > M is a 
continuous linear operator, F'(m) • u m — > F'(u) • w. Hence F'(m) • u > for 
any w g V. 
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V can be characterized in a more convenient way: 
Lemma 3. For any u E K, 

V = W := {v E X : Bv > on uj c }. 
We postpone the proof to the end of this section. 

V is a closed convex cone, with dual cone 

V ± = {7 G X' : (7, v) > VvE V}. 




Lemma 4. If f E P , then supp / C UJ C and f > 0. 

Proof. Given any y with support in y £ P and — y e P. Hence (/, y) = 
and therefore supp / C o> c . Now take y E Y positive. Then in particular 
y > on cj c , and y E P. By definition of this implies / > 0. 

We now use the following Lemma to characterize V in a different way. 

Lemma 5. Let Y be a Banach space, and P C Y a closed convex cone with 
dual cone P 1 - . Let X be a second Banach space, and A : X — > Y a bounded 
linear operator. Let K be the following cone in X : 

K = {uEX : AuE P}. 

Then the dual cone can be characterized by 



The proof of this Lemma can be found in [2] . An immediate consequence 
of Lemma 5 is 

Corollary 2. 



We now turn to the proof of the main theorem of this section. 

Proof (Proof of Theorem 2). We have seen that, since u is a minimizer, 
F'{u) E V ± and 

V = {v E X : B(v) > on uj c }, 

by Lemmas 2 and 3. By Corollary 2 there exists an / E P- 1 such that 
F'(u) = B T f, and by Lemma 4 supp / C oj c and / > 0. The conjugate 
operator B T is easily seen to be given by 



K 1 - = {A T g EX'-.gE P x }. 



V ={B T fEX':fEP ± }. 




(19) 
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for a smooth function <p E Y' , where <fi is implicitly extended by zero outside 
of the interval [0,T — 1]. We use the same notation for a general Radon 
measure f E Y' . 

Lastly, direct computation gives 

pT 

F'(u)-v = / [2a(u)uv +a(u)u' 2 v + b'(u)v], (20) 
Jo 

and hence we obtain the equation 

-2[a(u(x))u'(x)}' + a'(u(x))u'(x) 2 + b'(u(x)) = f f (21) 

Jx-l 

in the sense of distributions. 

We now turn to the statement of regularity. Since / G RM([0, T — 1]), 
the function g (defined in (5)) is uniformly bounded. Since all terms in (21) 
except the first are in L 1 (0,T), we have a(u)v! € W 1,1 , and the lower 
bound on a implies that u E VK 2,1 (0, T). Since W 2 ' 1 C W x -°° , the second 
term is now known to be in L°° , and again the lower bound on a is used to 
obtain u E W 2 '°°(§,T). This regularity of u implies that the distributional 
equation (21) is also satisfied almost everywhere. 

We still owe the reader the proof of Lemma 3. 

Proof (Proof of Lemma 3). V C W: Since B : X — > Y is continuous, 
W is closed, and therefore it suffices to show that V C W . Take any v E V 
and x E uj c . Then B(u + e n v)(x) > 0, and by definition of u c , Bu{x) = 0, 
implying that Bv(x) > 0. It follows that v E W. 

W C V: First consider w E W such that supp (Bw)- (the support of 
the negative part of Bw) is contained in u)% We claim that w E V, for which 
we have to show that there exists 

{£„} C M + , e n -» 0, such that B(u + e n w) > Vn E N. 

For x E uj c , Bu(x) = 0, and since Bw{x) > we have B(u + e n w)(x) > 
for any sequence {s n } C M + . For the complement uj c c1 note that since 
supp is compact and contained in the open set there exists S > 

such that Bu > S > on supp Hence, if e n < 5 WBwW^L , then 

Bu+e n Bw > on supp (Bw)- . Note that Bu > on [0,T-1], and Bw>0 
on (supp (Bw)-) c . This means that Bu + e n Bw > on Together with 
Bw > on u) c , this implies w E V. 

Finally, consider a general w E W. Fix a smooth function <p E X with 
(j> > on (0, T); note that B<p > c > 0. We approximate w by the function 
w £ := w + E(j>. We claim that supp (Bw e )- C w c c for sufficiently small e > 0. 
It then follows that w £ E V and w £ — > w, implying that w G V. 

To prove the claim, note that w G X C Hence is Lipschitz 
continuous, with Lipschitz constant 2||u;||ioo. Hence, for small enough e, 

Bw E (x) = Bw(x) + eBcj)(x) 

> Bw(y) + eBcf>(y) - 3| |w| \ L «, \x - y\. 
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Suppose Bw e (x) < and y G uj c . Then Bw(y) > 0, and by (22), 
-eB4>{y) > -3\\w\\ L oo\x - y\, 

or 

\x-y\> emy) 



3IM 



Therefore <i(supp (Bw)^, uj c ) > Ce for a suitable C > 0. Hence supp (Bw s ). 
C lo c c for small enough s, which proves the claim. 



6. Characterization of stationary points 

For this section we assume that the conditions of Theorem 2 are met. 

Lemma 6. Let u be a stationary point, and let g be defined as in (5). 

1. For all x G uj c , u(x) = u(x + 1) and u'(x) < u'{x + 1). 

2. If lu c contains an interval [xq,xi], then 

— u'(x) = u'(x + 1) for all x G (xo,Xi); 

— u"{x) = u"(x + 1) and g{x) = g(x + 1) for almost all x G (x ,Xi). 

This lemma imposes an interesting form of periodicity on the solution 
and the right-hand side g. Although the constraint is a non-local one, on 
an interval of contact of length L the solution actually only has the degrees 
of freedom of an interval of length one; the other values follow from this 
assertion. 

Proof. Since x G u) c , 

r x+l 

u = 0. 



/ 

Hence, since Bu(x) = f* +1 u G W 3 >°°, and Bu > 0, 

d f x+1 

= — / u = u(x + 1) — u(x), 
dx J x 

and 

d 2 r x+1 

0<-^J u = u'(x + l)-u'(x). 

If Bu = on [a;o,a;i], then the inequality above becomes an equality a.e. on 
the interior (a;o,a;i), implying that 

u'(x) = u'(x + 1) on (x , X\). 

The periodicity of u" and g now follow from (18). 
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Lemma 7. Let u be a stationary point, and assume that uj c contains an 
interval I. Then 

r x+l 

9 

J X 

is constant on Int /. 



/ 

J x 



Proof. By Lemma 6 u(x) = u(x + 1) for all x e /, and u"(x) = u"(x + 1) 
a.c. on I. In addition, u'(x) = u'(x + 1) for all x G Int/. Hence 

r-x + 1 

x i > I [-2a(u)u" - a'{u)u' 2 + b'(u)] (23) 

J X 

is constant on /. But by (18), (23) is equal to 

rx + l r-s rx+l 
/ / /= / ^ 

J X J 8 — 1 J X 

The following two lemmas and the theorem that follows are essential 
in determining the structure of the right-hand side g and therefore of the 
measure /. The main argument is the following. The function g has no reason 
to be monotonic; its derivative in x equals f(x) — f(x — 1), and although 
/ is a positive measure this difference may be of either sign. However, if 
for instance a left end point xq of lo c is flanked by a non-contact interval 
(xq — l,Xo), then the measure / is zero on that interval, and the function 
g is non-decreasing on (x 0l x a + 1). It is this argument, repeated from both 
sides, that allows us to determine completely the structure of the function 
g and the underlying measure /. 

Notation Let [xo,xi] C w c . Define 

p = xi — x (mod 1), (24) 

and 

P = min{n e N : n > Xi - x }. (25) 
Throughout the rest of this paper r is the translation operator defined 

by 

(tu){x) = u(x + 1). (26) 

Lemma 8. Let u be a stationary point, such that uo c contains an interval 
[xq,Xi]. Assume furthermore that 

supp /n (x - l,x ) = 0. (27) 

Then 

1. if xi -xoeN, then g does not decrease on each of the subintervals 



(x + i, x + i + 1), i = 0, 1, . . . , P; 
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2. if x\ — x £ N, then g does not decrease on each of the subintervals 
(x + i, x Q + i + 1), i = 0, 1, . . . , P - 1, 

nor does it on 

(x + P,xi + 1). 

Proof. On (x ,x + 1), 

g' = f-T-\f (2 = 7) />o, 

and therefore g is non-decreasing on (x ,x + 1)- By Lemma 6, g(x) = 
g(x + 1) for almost all x G (xo,xi). This implies that on each consecutive 
interval (xo + i, xq + i + 1), i = 1, . . . , P — 1, g docs not decrease. By the 
same reasoning, if x\ — x € N, then this also holds for (x + P, x + P+ 1) = 
(x\,Xi + 1). If not, then it holds for (x + P, x\ + 1). 

Remark 2. Let u be a stationary point. Define the mirror image v(x) = 
u(T - x), and h(x) = f(T -x-1). Then (v, h) solves 

f -2a(v)v" - a'(v)v' 2 + b'(v) = h, 
1 «(0) = v(T) = 1, 

and hence is also a stationary point. 

Applying Lemma 8 to (v, h) yields for (u, /): 

Lemma 9. Let u be a stationary point such that lo c contains an interval 
[xq,x\]. Assume furthermore that 

SUpp f(l (X! + 1,X! +2) = 0. 

1. if x\ — x € N, then g does not increase on each of the subintervals 

(x + i, x a + i + 1), i = 0, 1, ... , P; 

2. if Xi — xq N, then g does not increase on each of the subintervals 

(.x +p + i, x +p + i + 1), i = 0, 1, . . . ,P - 1, 
nor does it on 

(x ,x +p). 

To combine the previous two Lemmas, let 

X.^xo + i, i = 0,...,P, , 2g . 

y; =x +p + », * = o,...,p. 1 J 
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Theorem 3. Let u be a stationary point such that the contact set lo c con- 
tains an interval [xq,x\\. Suppose that 

supp/n{(a;o-l,a;o) U (an + l,an + 2)} = 0. 
Then there exists GgI such that 

1. if xi — x e N, then g = G on (x ,xi + 1), and 

p 

/l(x -l,x 1+ 2) = G^S{x-Xi). 

i=0 

2. if x\ — x £ N, then 

a(r)-i 91 0n \- Xi ^ * = 0,...,P, 

9[X) -\g 2 on[Y i ,X i+1 ],i = 0,...,P-l, ^> 



and 



p-i 



f\(x a -l,x 1+ 2) = ai§ ( X ~ X ') + b * S ( X ~ 

i=0 

where Oj = (G — -p)gi and bi = ^p^gi, and 

GP / PG 

9i ■■= Y 



.92 



GjP + l) P + l / (P + 1)G \ 

- p + i-p = e V G ' p J • (31) 



Proof. (1) h-ioGN. 

By Lemma 8, g does not decrease on the intervals (JQ, i = 
0, 1, . . . , P, and by Lemma 9 g does not increase on these intervals either. 
Hence g is constant on each interval. By Lemma 7 the constant is the same 
on each interval, i.e. that g = G on (x 0l xi + 1). This also implies that 
within the interval (xq — l,an + 2), / can only have support in the points 
x = Xo,X\, . . . , Xp = an, yielding the formula for / in the statement of 
the theorem. 

(2) x x - x N. 

Combining Lemma 8 and 9, we find that g is constant on each interval 
(Xi, Yi) and (Yi,X i+1 ), i = 0, 1, . . . , P - 1, and on (X P , Y P ). By Lemma 
6, g(x) = g(x + 1) for almost all x € (20, an), and hence g takes three 
values, 0, and g\ and g 2 (say) on (a:o,an + 1). We choose g = g\ on (Xi, Yi), 
i = 0,1,..., P, and g = g 2 on the intervals inbetween, (Yi, Xi+±), i = 
0, 1, . . . , P — 1; outside of the interval (xq, an), g vanishes. By Lemma 7, 

rx+l rx+p rx+1 

G= 9= gi+ g 2 = pgi + (1 - p)92- (32) 

J x J x J x+p 

Either g x = g 2 = G or g x < G < g 2 . The first case implies that g = G on 
[x ,xi + 1]. This implies that / does not only have support in X , X\, . . ., 
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Xp_i, but by reasoning for the mirror image (v, h) it also implies that / has 
support in Yq,Y\,... , Yp_i = x\. This is impossible. Hence g\ < G < g 2 . 
The support of / on (xo — 1, x\ + 2) is now seen be to limited to the set 
given in the statement of the Theorem. 

Thus we conclude that / is a sum of delta functions, but we still have 
to determine the weights a, and 6j. Since / = on (x — l,2?o)i we have 
gi = g(xo+) = f(xo)- Here we abuse notation, and write f(x) for the weight 
of the Dirac delta function at x. Now we have the following recurrence 
relations: 

f(X i ) + f(Y i )=g 2 , 
f(Y i ) + f(X i+1 )=g 1 , 

for i = 0, 1, . . . , P — 1. Solving this system we obtain 

f(Xi) = f(X ) - i(g 2 - gi) = .91 - i(g 2 - gi), 
f(Y i ) = (i + l)(g 2 -g 1 ). 

In addition, since x\ = Yp_i, f(x\) = P(g 2 — gi). On the other hand, 
gi = h(T — x\ — 1) = ,f{xi). This implies 

P + l 

.92 = —p—gi- 

To conclude, 

a i = f(X i ) = (G-±)g 1 , 

and 

bi = W) = % -±±g 1 . 

By Lemma 7, 

P9i + (1 -P)— p— 5i = G, 

which yields 

P = p(i--Wi. 

Solving for 171 now yields all required results. 

As we will see in the next Section, the contact set of u is connected 
in many important cases. Hence Theorem 3 allows us to give concise ex- 
pressions for g in cases that ui c is an interval of positive length (using the 
Heaviside function H): 

Corollary 3. If the contact set is an interval of positive length, then g 
equals the explicit function 

g(x;x ,x 1 ,G) = 

H(x — xo) — H(x — xi — 1) if x\ — xo e N, 

< 0i(JT(a;-a;o)-ir(a;-a;i-l))+ (33) 
+ (.92 - .9i) Ef=i [H(x - Xi) - H(x - Yi)] if Xl - x i N. 
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Here the coefficients gi, 2 are computed from Xq, X\, and G by (24), (25), 
(30). and (31). 

Figure 6 shows examples of both cases. For the remaining two cases of a 
stationary point that has a single or no contact point, g is immediately 
clear: with a single contact point, 

fm on [x ,x + 1], 
g(x) = < 

1 otherwise, 

for a suitable constant m > 0, while when there is no contact then obviously 
.9 = 0. 




Fig. 6. A generic picture of g(x;x ,T) for T — xo £ N (a) and for T — x G N 
(b). The light gray shaded areas represent the contributions of the individual delta 
functions of the corresponding Radon measures /. As an example, in (b) f consists 
of four Dirac deltas, all with mass G, at xo = Xq, xq + 1 = Xi, . . . , xo + 3 = X3. 
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7. The contact set is an interval 



In order to extract more information on the right-hand side g and the 
measure / than that given by Theorem 2 we study two cases. In the first 
case we assume that the operator given by the left-hand side in (18) satisfies 
a version of the classical comparison principle. In the second case we restrict 
ourselves to global minimizers. 

Definition 2. Let N be a (non)lincar operator on U. N is said to satisfy 
the comparison principle if for any [xq,xi] C [0, T], 



See e.g. [9] or [17] for a general exposition. Operators of the type considered 
here, i.e. 



may fail to satisfy the comparison principle for two reasons. First, the zero- 
order term b'(u) need not be increasing in u; for instance, the operator 
u i ► — u" — u does not satisfy the comparison principle on any interval of 
length 2ir or more. In a slightly more subtle manner, the prefactor a(u) of 
the second-order derivative may also invalidate the comparison principle; 
see e.g. [9, Section 10.3] for an example. 

We conjecture that the 'true' rod functions a and b given in (1) do not 
give rise to a comparison principle: b' is not monotonic, suggesting that on 
sufficiently large intervals the principle will fail. 

We first prove a lemma that will be used in both cases. 
Lemma 10. Let u be a stationary point such that x\,x 2 <G uj c . Assume that 




Nu := -2a(u)u" - a'(u)u' 2 + b'(u), 



(xi,x 2 ) n uj c = 0. 



Then 




(34) 



and for any m € (0, x 2 — Xi) 




and 




(35) 



Proof. Since (xi,x 2 ) n uo c = 0, 




u > 
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for all m G (0, x 2 — xi), which implies 

/•xi+m+l pxi+m pxi+1 />xi+m+l />xi+m 



/ u — / U = u+ u — / u 

<.xi+m+l 

= / u > 0. 

The other two assertions are handled similarly. 



Theorem 4. Let u be a stationary point, and assume that 

Nu -2a(u)u" - a'(u)u /2 + b'{u) (36) 

satisfies the comparison principle. Then uj c is connected. 

Proof. We proceed by contradiction. Since o> c is closed, non-connectedness 
implies the existence of x\, x 2 <G lo c such that (xi, x 2 ) H lo c = 0. 

Set v = u — tu. Then v(xi) = v(x2) — by Lemma 6, J" 2 v = by (34), 
and 

xi+m 

« < for all < m < x 2 - x\ (37) 

i 

by (35). Hence there exists anie (xi,X2) such that v(x) = 0. 

From (xi, x 2 ) fl u> c = it follows that supp / n (xi, x 2 ) = 0- Hence g = 
fx-i f 1S a decreasing function on (xi, x 2 ) and rg is an increasing function 
on this interval by previous arguments. Hence g — rg is a decreasing function 
on (xi, x 2 ). There arc three possibilities, each leading to a contradiction with 
the comparison principle. 

Case 1: g > rg on (xi, x 2 ). On (xi,x), 

Nu — g > rg — Ntu, 

u(Xi) = Tu(xi), 
u(x) — Tu(x). 

By the comparison principle, u > tu on (xi,x), i.e. v > 0. But this contra- 
dicts (37). 

Case 2: there exists an x such that g > rg on (xi,x) and g < rg on 
(x,x 2 ). If x > xi, the same argument applies. If x < x, we consider (x,x 2 ) 
instead, and apply the same argument. Now we conclude v < on (x,x 2 ). 
But observe that from f x v < by (37) and f X2 v = we have f X2 v > 0, 
which again implies a contradiction. 

Case 3: g < rg on (xi,x 2 ). Again we obtain a contradiction from con- 
sidering the interval (x,x 2 ). 



L 
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For the second case we limit ourselves to global minimizers. The results 
of this theorem do apply to the functions a and b given in (1). 

Theorem 5. Let u be a minimizer. Assume that a and b are of class C 1 
and that a is strictly positive. Then lu c is connected. 

Proof. As in the proof of Theorem 4 we assume that there exist x\,x 2 <G lo c 
with (xi, x 2 ) n uj c = to force a contradiction. Then 

supp fn(x u x 2 ) = 0, (38) 

and hence g is a decreasing function on (x\, x 2 ), and an increasing function 
on (x\ + l,x 2 + 1)- Now consider the following two new functions 

!u(x) on [0, xi], 
u(x + 1) on [xi,X2], 
u{x) on [x 2 ,T], 

and 

(u(x) on [0, xi + 1], 
u(x — 1) on [x\ + 1, x 2 + 1], 
u(x) on[x 2 + l,T]. 

Both are admissible, i.e. v,w e K: they are continuous by Lemma 6, imply- 
ing that v,w <G X, and the fact that Bv,Bw > follows from Lemma 10. 
In fact we need certain strict inequalities, which we derive after introducing 
some notation. 

The functions v and w are minimizers. To show this, write 
F{u\[ Xl , X2 ]) = [a(u)u' +b(u)]. 

J X\ 

Then since u is a minimizer, and since u and v only differ on [xi,x 2 ], 
F(u\ [xuX2] ) < F(v\ [xuX2] ) = F(w| [xi+liX2+ i]), 

and similarly 

^( U I[X1+1,X2+1]) < ^( W I[X1+1,X2+1]) = F ( U \[x U X 2 ])- 

This implies that 

^H[xi,x 2 ]) = - F ( u l[x 1 + l,x 2 + l]), 

and that F(u) — F(v) — F(w). Every minimizer is also a stationary point, 
and hence for v and w there exist positive Radon measures f v and /„, such 
that supp f v C lo c {v) and supp /„, C uj c (w). We also denote g v (x) — j x _ x f v 

and g w (x) = J x _ 1 f w - 
For any x € (xi,x 2 ), 

r x+l 

u > 0. (39) 



/ 

J X 
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Let first x 2 > x\ + 1. Then for any x E (x\ — l,x\), 

{•X+1 fX\ fX+1 



rX + L rXl rX + L 

/ V = V + V 

J x J X J x 1 

rxi r-x 

= U + 

J x J x-i 



x+2 

U 

+ 1 

rx\ rx + 1 

> / u+ u by Lemma 10 

J x Jxi 



-L 



x + 1 

U > 0. 



For any x <G (xi, x 2 — 1) the same is true: 

f-x + l r x+2 



rx+L rx+i 

/ ^ / 

Jx J x+1 



U > 0, 



since x+1 < x 2 , which allows us to use (39). Now let x 2 < x\ + 1. Then for 
any x £ (xi — 1, x 2 — 1), we can repeat the first argument above to conclude 



I 



x+1 

v>0. 



Combining these statements we find 

rx+l 

I v > for all x e (x\ — l,x 2 — 1), 

J X 

which implies uj c {v)f\{xi — 1, x 2 — 1) = 0. Hence supp f v f\{xi — l,x 2 — 1) = 0. 
But since u and v conincide on [0, £1], we have (/t;|[o.a;i] = 9u\[o,xi]^ so that 
supp / M n(a;i — 1, x 2 — 1) = 0. Combined with (38), this implies that g u \[xi,x 2 ] 
is constant. By symmetry the same is true for gu\[ Xl +i.x 2 +i]- Note that if 
x 2 > x\ + 1, then the overlap implies that the two constants are the same; 
for the other case we now prove this. 

Define z = u — tu; the function z solves the equation 

- 2a(u)z" = g u - rg u + {a {u)u' 2 - a'(Tu)(™)' 2 } 

- {b'(u) - b'(Tu)} + {2a(u) - a(Tu)}(Tu)" (40) 

on the interval (xi,x 2 ). Of the right-hand side, we have seen above that 
the term g u — rg u is constant on (xi, x 2 ); let us suppose it non-zero for the 
purpose of contradiction. The function z is of class C , and both z and z' 
vanish at x = Xi :2 . Therefore the assumed regularity on a and b implies 
that the expressions between braces are continuous on [x i,x 2 ] and zero at 
x = xi, 2 . The sign of the right-hand side of (40) is therefore determined 
by g u — rg u , and most importantly, is the same at both ends x\ and x 2 ; 
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therefore the sign of z, at x = Xi + and x = x 2 — , is also the same. This 
contradicts the following consequence of Lemma 10: 



«ltm rx 2 

I z < and / z > for all < m < x 2 — X\. 

J Xi J X2—m 

This leaves g u = rg u on (0:1,0:2). But then, by uniqueness of the initial- 
value problem, u(x) = u(x + 1) for all x G [0:1,0:2], and [a; 1,0:2] C co c , 
contrary to our assumption that (0:1, 0:2) n uj c = 0. 



8. Symmetry 

In the introduction we raised the question whether the stationary points 
or minimizcrs inherit the symmetry of the formulation, or to put it differ- 
ently, whether non-symmetric solutions exist. 

For the discussion of this question it is useful to introduce an equivalent 
formulation of the Euler-Lagrange equation (18) similar to the Hamiltonian- 
systems formulation used in the proof of Theorem 1 . For the length of this 
section we assume that Theorem 3 applies and therefore that there is a 
single contact interval [20,211]. 

By multiplying (18) with v! and integrating one finds that the function 
H , defined by 

H := -a{u)u' 2 + b(u) - gu, (41) 

is piecewise constant, and that H and g jump at the same values of x. The 
function g takes three values on [0, T], these being the values g\ and g 2 
introduced in Theorem 3, and the value go = outside of the extended 
contact interval [0:0, a; 1 + 1]. (Note that g\ and gi may be equal). We claim 
that H also takes three values, Ho, H\, and H2, and that these values 
correspond to those of g, i.e. that the pair (g, H) takes three values (0, Ho), 
{gi,Hi), and (g 2 ,H 2 ) (although it may happen that (#i,ifi) = (g 2 ,H 2 )). 

To prove this claim, first consider the case of p > 0, where p is defined 
as in (24). Then 

u\( Xo ,x a +p) = m|(x +1,x +1+p) an d g\(x ,x +p) = <7|(x +l,x +l+p) 

by Lemma 6 and (29). Therefore H is the same on these two intervals. 
Repeating this argument for all subintervals of [xo, x\ + 1] of the form (xo + 
k,x + k+p) and (x + k +p, x + k + 1) we find that H takes two values on 
the interval [xo, x\ + 1], Hi and H 2 , and that these coincide with the values 
gx and g 2 of g. 

When p = 0, a similar argument yields that H takes only one value on 
[x ,xi + 1] (as does g). 

A consequence of this characterization of H is the following lemma: 
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Lemma 11. Under the conditions and notation of Theorem 3, 

u(x ) = u(x + p) = u(x + 1) = u(x + 1 + p) = ■ ■ ■ = u(xi + 1). 

Proof. When p = the statement follows from Lemma 6. For p > 0, note 
that at any of the interior jump points, i.e. at all jump points except xq 
and x\ + 1, we have [H] — ~[g]u where [H] = ±(i?2 — Hi) and [g] = 
±(<?2 — Regardless of the sign this equation has only one solution u. For 
the remaining two points x and x\ + 1 the result follows from Lemma 6. 

We still need to show that the value of H is the same on both sides of 
the extended contact interval [2:0,2:1 + 1], so that we can define the value 
H unambiguously. If one of the ends of this interval equals or T there is 
nothing to prove; we therefore assume that min{x , T — x\ — 1} > d > 0. 
Now multiply (18) with the function 



Since H is constant on (0, d) and on (T — d, T) the two constant values are 
equal; we then define H to be this value. 

We now turn to the implications of this characterization of solutions 
(u,g) and the associated pseudo-Hamiltonian function H. 

Theorem 6. Let u be a stationary point with a single contact interval [xq, x±] . 
Let p be given as in (24), and define the set of jump points J = {x n ,x + 
p,x + l,x a + l+ p,...,x 1 + l}. 

1. There exists a£l such that at any x e J, u'(x) = ±a. 
Now assume that b is non- decreasing on [l,oo). 

2. If the operator N given in (36) satisfies the comparison principle, then 
u is symmetric on [0, T]. 

3. If u is a minimizer with u'(x ) = —u'(x\ + 1), then u is symmetric on 



v(x) 



>'(*) 
< u'(x) 



< x < d 
d<x<T-d 
T - d < x <T, 



{^u'(x) 



and integrate to find 




[0,T]. 



Proof. For the first part write 



u 



b{u) — gu — H 



a(u) 



and note that by the proof of Lemma 11 the sum gu + H is continuous. 
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For the second part, note that by Lemma 11m has the same value on 
each end of the interval [xo, %o + p] (if p > 0) or [xo, x n + 1] (if p = 0). By 
the uniqueness that follows from the comparison principle the function u is 
symmetric on this interval. By repeating this argument over all subintervals 
of [xq, x\ + 1] we find that u is symmetric on [xq,x\ + 1]. 

The functions u\{t) := u(x — t) and u 2 (t) := u{x\ + 1 + t), therefore, 
have the same zeroth and first derivatives at t = 0; they satisfy the same 
equation (41) (note that H is symmetric on [x , x\ + 1]); therefore the two 
functions are equal as long as they both exist. This implies that lack of 
symmetry must stem from a difference in domain of definition of Mi and m 2 



We claim that neither u\ nor u 2 has an interior maximum. Assuming 
this claim, the assertion of the theorem follows since the monotonicity of 
mi, 2 then implies that the boundary condition 1*1,2(2) = 1 has at most one 
solution t. 

Now assume that u\ has a maximum at t\ > 0. The function Mi is 
solution of the Hamiltonian system (41), where H and g are constant for 
t > 0. As in the proof of Theorem 1, therefore iti(ti) > 1. Choose a bounded 
interval / C [0, 00) such that u > 1 on Int/ and u(dl) = 1. 

The reduced functional F(v) = Jj[a(v)v' +b(v)] has a global minimizcr 
v in the class of functions v satisfying v{dl) — 1. From studying the pertur- 
bation v 1 ► min{u, 1} and using the monotonicity of b it follows that v < 1 
on /. By the comparison principle this is the only stationary point of F, a 
conclusion that contradicts the fact that u\ is a different stationary point. 

For the third part, first note that the support of the continuous function 
x 1 ► J^ +1 u is [xq, Xi]; therefore 

for every e > there exists S > such that any perturbation v with 
d(supp v, {x n7 xi + 1]) > e is admissible provided < ^. 

We will use this below. 

The assumption on the derivatives places us in the same position as 
above: the functions u\(t) := m(x — t) and u 2 (t) :— u(x\ + l + t) are equal 
as long as they both exist. Again we will show that neither may have an 
interior maximum, but by a different argument. 

Assume that u\ has a maximum. By defining t\ = x the boundary 
condition on u takes the form Ui(ti) = 1. Pick 



and define the set S = {t e [0,ii] : u\(t) > 0); we can assume that for 
e = inf S > we have max{ui(£) : < t < t{\ — (3 < S for the associated 5 
given above. 

Now define v(t) = mm{/3, u\(t)}. The function v is admissible by con- 
struction; it differs from Mi only on the set S, and therefore the difference 
in energy is given by (with a slight abuse of notation) 



for t > 0. 



max{l,Mi(0)} < [3 < max{Mi(£) : < t < t{\ 
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This contradicts the assumption of minimality. 

The conditions of Theorem 6 are quite sharp. We demonstrate this with 
two examples. 

Example 1: b is decreasing on [1, oo). It is relatively straightforward 
to construct a non-symmetric stationary point by choosing an appropriate 
function 6 that is decreasing on [l,oo), thus showing that part 2 of Theo- 
rem 6 is sharp. 

Take a symmetric stationary point u for which u < 1 on [0, T], u'{T) > 0, 
and for which the contact set is bounded away from x = T (see the next 
section for examples). Close to x = T, the function u satisfies 

,2 b(u) - H 
u = — — 

a{u) 

for some H S M., and since u'(T) > 0, 6(1) > H . Now change b{u) for u > 1 
such as to have (for instance) 6(2) = H, and continue the solution u past 
x = T. By construction u(T+f) = 2, for some f > 0, and u'(T+f) = 0; by 
symmetry then u(T + 2T) = 1. The new function u defined on the domain 
[0, t + 2T] is a non-symmetric stationary point (Figure 7). 




T T + 2T 

Fig. 7. A non-symmetric stationary point can be constructed by defining b(u) 
appropriately for u > 1. 



Example 2: equal (non-opposite) derivatives on du) c . For cer- 
tain functions 6 and domains [0, T] global minimization favours breaking of 
symmetry. We demonstrate this for the functional 

F(u) = J [u' 2 + a(l-u 2 ) 2 ], 

where a will be chosen appropriately. We consider the functional F on 
functions u : [0, 1] — > M with boundary conditions u(0) = u(l) = 0; although 
this is slightly different from the setup in the rest of the paper, it simplifies 
the argument, and the extension to a more general situation is intuitively 
clear. 
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We will show that 



inf |f(u): y* u >o| < inf <|V(u) : J u > and u is symmetric | . (42) 
The infimum on the right-hand side is bounded from below, 

F{u) > (1 - ac/2) j u' 2 + a, 
by the Poincare inequality 

u 2 < cd 2 / u' 2 for all u with u(0) = and I u = 0. (43) 
Jo Jo J 

The function v(x) = a + cos(6(l — x/d)) is optimal in this inequality, where 
a ~ 0.22 and 6 ~ 4.49 are determined by the boundary condition v(0) = 
and the integral condition J v = 0. The Poincare constant equals c ~ 0.0495. 
Note that for symmetric functions u we may take d = 1/2. 

At the function w(x) — sin27ra; the functional F has the value F(w) — 
2tt 2 + 3a/8. For all a G (16tt 2 /5, 2/c] ~ (31.6,40.3] therefore 

= 2-7T 2 + 3a/8 < a < inf |f(m) : J u > and u is symmetric| , 

which demonstrates (42). 

The reason for this preference for asymmetry can be recognized in the 
constant in the Poincare inequality (43) (sec Figure 8). For symmetric func- 
tions the relevant class is {u : [0, 1/2] — > M. : u(0) = J u = 0}, and for more 
general functions {u : [0, 1] — ► K : u(0) = u(l) = j u = 0}. For this latter 
class the Poincare coefficient is achieved by the function w above with the 
value c = 1/4tt 2 ~ 0.0253, which is larger than c(l/2) 2 = 0.0124. 





Fig. 8. Under symmetry conditions the effective domain, the domain on which 
J u — 0, is half the actual domain size. Equivalently, more (costly) oscillations are 
necessary. 
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9. Numerical simulations 

In this section we describe in detail our numerical simulations of sta- 
tionary points of F under constraint, i.e. of solutions of 



2a{u)u" - a'{u)u' 2 + b'(u) = f /, 

Jx-1 



(44) 



w(0) = u(T) = 1, (45) 
supp / C u> c , (46) 
/ a positive Radon measure, (47) 

/ u>0 Vxe[0,T-l]. (48) 

J X 

We concentrate on the case in which the solution is symmetric and the 
contact set is non-empty, and we use the fact that the right-hand side in the 
differential equation can be characterised explicitly (see (33)). We further 
simplify by replacing the inequality (48) by the condition that the function 
x i — ► f* +1 u has a second-degree zero at x — x , leading to the new system 
in the unknowns (u, x , G) 

- 2a(u)u" - a'{u)u' 2 + b'(u) = g(x; x , T-x Q -l,G), (49) 
u(0) = u{T) = 1, (50) 
u{x ) = u(x + 1), (51) 

r-xo + l 

/ u = 0. (52) 

J Xq 

For brevity we shall write g(x; xq,T 7 G) for g(x; x , T — x — 1, G). 

Lemma 12. Assume that the operator on the left-hand side of (4-9) satisfies 
the comparison principle. Then any solution of problem (44~4^) with non- 
empty contact set is also a solution of (49-52); vice versa, any solution 
of (49-52) is also a solution of (44~48)- 

Proof. Since the implication (44-48) => (49-52) follows by construction, 
it suffices to show the opposite implication; in fact, since an admissible 
measure / can be constructed from any g(x; xq, T, G), it is sufficient to 
show that solutions of (49-52) satisfy 



L 



x+1 

u>0 Vxe[0,T-l]. 



We show slightly more, namely that 

(■x+1 



J X 



U = Vie [XQ,T- Xq - 1] 
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and that 

rx+l 

/ u>0 Vx £ [x ,T- x - 1]. 

J X 

The function u is symmetric by Theorem 6. Since u(xq) — u(xq + 1), 

u(xo) — u(x + 1) = u(T — x ) = u(T — x — 1) =: u. 

Set u\{x) = u(xq + x), and u 2 {x) ~ u(xo + x + 1) for all x € [0, T— 2xo — 1]. 
By construction, g(x;xo,T,G) — g(x+ 1; xo, T, G) for all x G [x , T— x — 1]. 
Hence, if we set h(x) — g(x + x ; x , T, G), for all x e [0, T — 2x — 1], then 
Mi and u 2 both satisfy 

-2a(v)v" - a'(v)v' 2 + b'(v) = h, 
v(0) = v(T-2x Q -l) = u, 

By uniqueness, mi = u 2 on [0, T — 2x — 1]. In terms of u this means 
u(x) = u(x + 1) for all x G [x , T — x — 1]. But that implies that 

rx+l 

/ u = Vx € [xo, T — xo — 1]. 

It remains to be shown that 

/ u>0 Vx ^ [x ,r-x - 1]. (53) 

By symmetry we only show this for x < x . Let u p and g p be the 1-periodic 
extrapolation of u|[ Xo>Xo+ i] and g\[ XOtXo+ iy, note that f* + u p = for every 
x. For x < xo, 

-2a(u)u"-a'(u)u' 2 +6'(u) = < g p = -2a(u p )u' p -a! \u p )u' p 2 +b' (u p ), (54) 

implying that u > u p for x = x — and therefore also (53) for x = x — . 
If u and u p intersect again at some x < Xo, then the comparison principle 
and (54) imply that u < u p on [x, xo], in contradiction with the previous 
statement. This concludes the proof. 

We discuss two different ways of calculating solutions of the problem (49- 
52). 



9.1. Continuation 

We implemented a strategy of continuation of solutions, using the con- 
tinuation package AUTO [7], and we chose the simple case 



<*(«) = ^. b{u)= l -(u+lf. 



(55) 



36 



Van der Heijden, Peletier, and Planque 



To implement system (49-52) in AUTO, we divide [0, T] into three subdo- 
mains, [0,a:o], [a^o^o + 1] and [xq + 1,T] and specify the equations 



-u"{xi) +u 1 (x 1 ) + 1 = g(xi;x ,T,G) 

x[ = 1 

-U2(x 2 ) + u 2 (x 2 ) + 1 = g{x 2 ;x ,T, G) 

4 = 1 

-u'z(xz) + u 3 (x 3 ) + 1 = g(x 3l x ,T, G), 

4 = i 



on [0,xo], (56) 
on [x ,x + 1], (57) 
on[.T + l,T], (58) 



with boundary conditions 



ui(0) = 1, 
U!(x ) = u 2 (x ), u[(x ) = u' 2 (x ), 
u 2 (x ) = u 3 (x ), u 2 (x ) = u' 3 (x ), (59) 
us(T) = 1, 
u 2 (x ) = u 2 (x + 1), 
a;i(0) = 0, x 2 (x ) = x , x 3 (x + 1) = x + 1. 

and integral condition 



i X 



u 2 = 0. (60) 

I Xq 

Note that in (56)-(58) we have added trivial equations in order to solve for 
the Xi variables, which are required in the evaluation of g(x; xq, G, T). 

There are still some technicalities that have to be overcome: AUTO is not 
well-equipped to handle systems with a discontinuous right-hand side, such 
as the function g(x; x , G, T) that is supplied here. We remedy this by using 
a low-order method for all simulations, and we smooth the function g given 
in (33) by substituting arctans for Heaviside functions: 



g(x;x ,T, G) = 

^ ( arctan(,4(a; - x )) - arctan(A(x -T-x ))) if T - 2x Q <£ N, 
" Ef=i [arctan(A(x - X,)) - arctan(A( 2 ; - Y$)] , 

i ( arctan(A(x - x )) - arctan(A(x - T - x ))) if T - 2x e N, 

where Xi and Yi are as in (28). In the limit A — > oo, g(x: x , T, G) converges 
pointwise to g(x; Xo,T, G). 

There are nine differential equations with ten boundary conditions and 
one integral condition. This means that we expect to specify three free 
parameters to obtain a one-parameter curve of solutions. These are T, x , 
and an additional parameter (3. It worked well to choose the freedom in f3 
in modulating the values of g\. 2 \ 



9i =9i+ P and g 2 = g 2 + (3. 
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One may prove a priori that j3 = by remarking that 

r-xo + l rxo + 1 

/ 9= [-u" + u+l] =1, 

J x J x 

and using (32) to find 

rxo + l 

1= / g = p(gi+(3) + (l-p)(g2 + P) = l + 0. 

J x 

We have found no other role for (3 than to accommodate for small numerical 
inaccuracies due to the discontinuous right-hand side. In all simulations 
/3~ 1(T 4 . 

We have validated the code by comparing solutions from AUTO with 
explicit solutions. An example is given in Figure 9. 



1.2 



0.2 







-0.2 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



X 

Fig. 9. A comparison of a solution of system (56)-(60) produced with AUTO (o 
symbols) to an explicit solution, for a generic value of T (here scaled to 1): T = 
4.91635. In this simulation A = 1000. 

As we have seen in the discussion at the beginning of this Section, as T 
becomes larger the minimizer u has to have a contact point, and for large 
enough values even a full interval of contact. The point x , the leftmost point 
of contact, is determined as part of the solution; one may wonder how this 
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point depends on T. For operators N that satisfy the comparison principle, 
it is straightforward to prove that xo remains bounded for all T. Moreover, 
for the operator under consideration here, as T — > oo, x a — > log(2 + v / 3). 
These two phenomena are illustrated in Figure 10. 



1.33 



1.32 - 




1.26 1 1 1 1 1 1 1 

4 6 8 10 12 14 16 18 

Fig. 10. Behaviour of xo as a function of domain size T for system (56)-(60) 
computed with AUTO. As T grows, xo remains bounded and converges to log(2+\/3) 
(horizontal line). 



Since \g\ — g^\ — > as P (and therefore T) increases, g becomes con- 
stant in the limit of large T. By the comparison principle, u does the same, 
implying that F{u)/T — > 1. The start of the convergence to 1 is shown in 
Figure 11. 

9.2. Directly solving the boundary-value problem 

Computing solutions of the rod equations — rather than the simpler prob- 
lem (55) — using AUTO has proved difficult, for reasons that we do not un- 
derstand well. Instead, a boundary-value problem solver from Matlab was 
used to create Figure 3. Set 

^ 2u" 5uu' 2 3m a 

~47r 2 (l + u 2 )l + 4tt 2 (1 + u 2 )5 ~ r 2 (l + u 2 )l + (l + u 2 )f ' 
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1.125 




1.085 



Fig. 11. Behaviour of F(u) — T/2 as a function of domain size T for system 
(56)-(60) computed with AUTO. As T grows, F(u) - T/2 oscillates towards 1, the 
energy of u = on a unit length interval. 



To find a solution of 



Lu = g(x;x ,T,G), u(x ) = u(x + 1), / u = 0, 

Jxq 



for a generic value of T (large enough) we construct a two-parameter shoot- 
ing problem. Fix G and Xq and consider the boundary- value problem 





= 


on 


[0,x ], 


Lu 2 


= .9i 


on 


[x ,x +p], 


Lu 3 


= .92 


on 


[x +p,x + l], 


Lu u 


= .91 


on 


[xq + l,x + 1 +p], 


Lu 5 


= 


on 


[x + l +p,f], 



(61) 
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with boundary conditions 

«i(0) = 1, (62) 

U!(x ) = u 2 {x ), u[(x ) = u' 2 (x ), (63) 

u 2 (x Q + p) = u 3 (x Q + p), u' 2 (x + p) = u' 3 (x + p) , (64) 

"3(2:0 + 1) = u A (x + 1), u' 3 (x + 1) = u' 4 (x + 1), (65) 
ii 4 (x + 1 + p) = u 5 (x Q + l+ p), u' 4 (x + 1 +p) = u' 5 (x + l+ p), (66) 

u 5 (f) = 1. (67) 



Here, as before, p = T-2x -l (mod 1), P = min{n e N : n > T-2a; -l}, 
and 

GP _ G(P + l) 

91 ~ P+l-p' 92 ~ P+l-p' 

by Theorem 3. Note that this is not exactly the same problem as (49-50), 
since the periodic section has been reduced from P periods to a single period, 
and the solution is defined correspondingly on a smaller domain of length 

f = 2x +p+l. 

This allows us to use the decomposition in five subdomains for any T, which 
facilitates computation. This is illustrated in Figure 12. 



u 








T 


\ T 



h 

T - 2x -p-l 



Fig. 12. Schematic picture of the idea behind T — 2xo + p + 1. Since u (solid 
black line) is periodic between xo and T — xo, we can cut out an interval of length 
T — 2xo — p — 1 and find the corresponding solution on [0, T]. 

We now vary x and G to find solutions of system (61)-(67) that satisfy 

rx +p rxo+l 
U2{x ) = u 3 (x + 1), / U 2 + U 3 = 0, 

J X J XQ+p 

using a standard Matlab boundary-value problem solver, bvp4c. An example 
solution is drawn in Figure 3 in which we have used a = l/2n, r = 1. Note 
that all analysis in this paper assumes zero rod thickness; in Figure 3 the 
rod has been artificially fattened for better viewing. 
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